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Abstract. We use new algorithms, based on the finite lattice method of series 
expansion, to extend the enumeration of self-avoiding walks and polygons on the 
triangular lattice to length 40 and 60, respectively. For self-avoiding walks to length 40 
we also calculate series for the metric properties of mean-square end-to-end distance, 
mean-square radius of gyration and the mean-square distance of a monomer from the 
end points. For self-avoiding polygons to length 58 we calculate series for the mean- 
square radius of gyration and the first 10 moments of the area. Analysis of the series 
yields accurate estimates for the connective constant of triangular self-avoiding walks, 
fi = 4.150797226(26), and confirms to a high degree of accuracy several theoretical 
predictions for universal critical exponents and amplitude combinations. 
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1. Introduction 

Self-avoiding walks (SAWs) and polygons (SAPs) on regular lattices are combinatorial 
problems of tremendous inherent interest as well as serving as simple models of polymers 
and vesicles |25 | ITo ] IT6] . They are fundamental problems in lattice statistical mechanics. 
An n-step self-avoiding walk a; is a sequence of distinct vertices uq, wi, ■ ■ ■ ,&n such that 
each vertex is a nearest neighbour of it predecessor. SAWs are considered distinct up to 
translations of the starting point Uq. We shall use the symbol fl n to mean the set of all 
SAWs of length n. A self-avoiding polygon of length nisan - 1-st ep SAW such that ujq 
and u n -i are nearest neighbours and a closed loop can be formed by inserting a single 
additional step. One is interested in the number of SAWs and SAPs of length n, various 
metric properties such as the radius of gyration, and for SAPs one can also ask about 
the area enclosed by the polygon. In this paper we study the following properties: 

(a) the number of n-step self- avoiding walks c n ; 

(b) the number of n-step self-avoiding polygons p n ; 

(c) the mean-square end-to-end distance of n-step SAWs {R 2 ) n ] 

(d) the mean-square radius of gyration of n-step SAWs (R 2 ) n ', 

(e) the mean-square distance of a monomer from the end points of n-step SAWs {R, 

(f) the mean-square radius of gyration of n-step SAPs (R 2 ) n ', and 

(g) the k th moment of the area of n-step SAPs (a k ) n . 

The metric properties for SAWs are defined by, 



(Rl)n = — -^n) 2 , 

i r i n 

n n„ L v ; m=o 

i r i n 



u 3 f] 



with a similar definition for the radius of gyration of SAPs. 

It is generally believed that the quantities listed above has the asymptotic forms 
as n — » oo: 



c n = A f i n n'y- 1 [l + o(l)], 


(1.1a) 


p n = B^n a - 3 [l + o(l)], 


(1.16) 


(Rl) n = Cn^[l + o(l)l 


(1.1c) 


(Rl) n = Dn^[l + o(l)}, 


(l.ld) 


(R 2 m } n = En 2 »[l + o(l)), 


(Lie) 


(R\ = Fn 2u [l + o{l)], 


(1.1)) 


(a k ) n = G {k) n 2vk [l + o(l)]. 


(1-1^) 
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The critical exponents are believed to be universal in that they only depend on 
the dimension of the underlying lattice. The connective constant \x and the critical 
amplitudes A—G^ vary from lattice to lattice. In two dimensions the critical exponents 
7 = 43/32, a = 1/2 and v = 3/4 have been predicted exactly, though non-rigorously, 
using Coulomb-gas arguments [27^ l2"?j. 

While the amplitudes are non-universal, there are many universal amplitude 
combinations. Any ratio of the metric SAW amplitudes, e.g. D/C and E/C, is expected 
to be universal |Hj • Of particular interest is the linear combination [Bl E] (which we shall 
call the CSCPS relation) 

where y t = 1/v and yh = 1 + r y/{2v). In two dimensions Cardy and Saleur (as 
corrected by Caracciolo, Pelissetto and Sokal 0) have predicted, using conformal field 
theory, that H = 0. This conclusion has been confirmed by previous high-precision 
Monte Carlo work [2] as well as by series extrapolations [14] . 

Privman and Redner [2H| proved that the combination BC/aao is universal, Cardy 
and Guttmann [I] proved that BF = -^^aao, and Cardy and Mussardo [H| proved that 
C/F is universal and gave the first theoretical estimate of the value C / F 13.70. a 
is an integer constant such that p n is non-zero when n is divisible by a. So a = 1 for 
the triangular lattice and 2 for the square and honeycomb lattices, ao is the area per 
lattice site and ao = 1 for the square lattice, ao = 3^/3/4 for the honeycomb lattice, and 
a = a/3/2 for the triangular lattice. 

Richard, Guttmann and Jensen [SUj conjectured the exact form of the critical 
scaling function for self-avoiding polygons and consequently showed that the amplitude 
combinations G^B^ 1 are universal and predicted their exact values. The exact value 
for G^ = j- had previously been predicted by Cardy [3]. The validity of this conjecture 
was recently confirmed numerically to a high degree of accuracy using exact enumeration 
data for SAPs on the square, honeycomb, and triangular lattices |3*T| . 

The asymptotic form (jl.laj) only explicitly gives the leading contribution. In general 
one would expect corrections to scaling so, e.g, 



A[i n n 



a\ a 2 b bi b 2 



n n 2 n Al n Al+1 n Al 



1.3) 



In addition to "analytic" corrections to scaling of the form ak/n h , there are "non- 
analytic" corrections to scaling of the form bk/n Al+k , where the correction-to-scaling 
exponent Ai isn't an integer. In fact one would expect a whole sequence of correction- 
to-scaling exponents A x < A 2 . . ., which are both universal and also independent of the 
observable, that is, the same for c n , p n , and so on. In a recent paper pQ we study the 
amplitudes and the correction-to-scaling exponents for two-dimensional SAWs, using 
a combination of series-extrapolation and Monte Carlo methods. We enumerated all 
self-avoiding walks up to 59 steps on the square lattice, and up to 40 steps on the 
triangular lattice, measuring the metric properties mentioned above, and then carried 
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out a detailed and careful analysis of the data in order to accurately estimate the 
amplitudes and correction-to-scaling exponent. The analysis provides firm numerical 
evidence that Ai = 3/2 as predicted by Nienhuis [2l)| I27j. 

In this paper we give a detailed account of the algorithm used to calculate the 
triangular lattice series analysed in PIE], perform some further analysis of the series 
and confirm to great accuracy the predicted exact values of the critical exponents, then 
we briefly summarise the results of the analysis from QjJ |M] and finally study other 
amplitude combinations. 

2. Enumeration of self-avoiding walks and polygons 

The use of transfer-matrix methods for the enumeration of lattice objects has its origin in 
the pioneering work of Enting 1.0] who enumerated square lattice self-avoiding polygons 
using the finite lattice method. The basic idea of the finite lattice method is to calculate 
partial generating functions for various properties of a given model on finite pieces, say 
W x L rectangles of the square lattice, and then reconstruct a series expansion for the 
infinite lattice limit by combining the results from the finite pieces. The generating 
function for any finite piece is calculated using transfer matrix (TM) techniques. The 
algorithm we use to enumerate SAWs and SAPs on the triangular lattice builds on this 
approach and more specifically our algorithm is based in large part on the one devised 
by Enting and Guttmann ^1] for the enumeration of SAPs on the triangular lattice 
with the generalisation to SAWs following the work of Conway, Enting and Guttmann 
[Zj and using further recent enhancements and parallelisation as described in [TOJ ED] • 

2.1. Basic transfer matrix algorithm 

In this section we give a detailed description of the SAW algorithm and then briefly 
outline the changes required to enumerate SAPs. 

We implement the triangular lattice as a square lattice with additional edges 
connecting the top- left and bottom-right vertices of each unit cell (see figHJ). We use 
W x L rectangles as our finite lattices. The most efficient implementation of the TM 
algorithm generally involves bisecting the finite lattice with a boundary (this is just a 
line in the case of rectangles) and moving the boundary in such a way as to build up the 
lattice cell by cell. The sum over all contributing graphs is calculated as the boundary is 
moved through the lattice. Due to the symmetry of the triangular lattice we need only 
consider rectangles with L > W. SAWs in a given rectangle are enumerated by moving 
the intersection so as to add one vertex at a time, as shown in figure H In most cases it 
is most efficient to let the boundary line cut through the edges of the lattice. However, 
on the triangular lattice it is more efficient to let the boundary line cut through the 
vertices [H]. Essentially this variation leads to only half as many intersected vertices 
(as opposed to edges) along the boundary line. For each configuration of occupied or 
empty vertices along the intersection we maintain a generating function for partial walks 
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Figure 1. A snapshot of the boundary line (dashed line) during the transfer matrix 
calculation on the triangular lattice. SAWs are enumerated by successive moves of 
the kink in the boundary line so that one vertex (shaded) at a time is added to the 
rectangle. To the left of the boundary line we have drawn an example of a partially 
completed SAW. 

cutting the intersection in that particular pattern. If we draw a SAW and then cut it 
by a line we observe that the partial SAW to the left of this line consists of a number 
of loops connecting two vertices (we shall refer to these vertices as loop-ends) in the 
intersection, and pieces which are connected to only one vertex (we call these free ends). 
The other end of the free piece is an end point of the SAW so there are at most two free 
ends. In addition it is possible that the SAW touches a vertex (that is the SAW comes in 
along one edge and exits along another edge both without crossing the boundary line). 
All these cases are illustrated in figure ^ In applying the transfer matrix technique to 
the enumeration of SAWs we regard them as sets of edges on the finite lattice with the 
properties: 

(1) A weight u is associated with an occupied edge. 

(2) All vertices are of degree 0, 1 or 2. 

(3) There are at most two vertices of degree 1 and the final graph has exactly two 
vertices of degree 1 (the end points of the SAW). 

(4) Apart from isolated sites, the final graph has a single connected component. 

(5) Each graph must span the rectangle from left to right and from bottom to top. 

We are not allowed to form closed loops, so two loop-ends can only be joined if they 
belong to different loops. To exclude loops which close on themselves we need to label 
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Figure 2. The seven possible outputs from a single iteration of the TM algorithm. 
Depending on the states of the three vertices a, b, and d in the input some of the 
outputs cannot occur. 



the occupied vertices in such a way that we can easily determine whether or not two 
loop-ends belong to the same loop. The most obvious choice would be to give each loop 
a unique label. However, on two-dimensional lattices there is a more compact scheme 
relying on the fact that two loops can never intertwine. Each end of a loop is assigned 
one of two labels depending on whether it is the lower end or the upper end of a loop. 
Each configuration along the boundary line can thus be represented by a set of states 
{(Tj}, where 



empty vertex, 

1 vertex is a lower loop-end, 

2 vertex is an upper loop-end, 

3 touched (degree 2) vertex, 

4 vertex is a free end. 



'XI) 



If we read from the bottom to the top, the configuration along the intersection of the 
partial SAW in figured is {310341022}. 



2.1.1. Updating rules In figur^21 we have illustrated what can happen locally as the 
boundary line is moved. Before the move, the boundary line intersects the vertices a, b 
and d and after the move the vertices a, c and d are intersected by the boundary line. 
We shall refer to the boundary line configuration prior to a move as the 'source' and 
after the move as the 'target'. In a basic iteration step we can insert bonds along the 
edges emanating from vertex b. Since vertex b can't have degree greater than 2 we can 
at most insert two new bonds. However, depending on the states of vertices a and d in 
the source, some of the edge configuration in figure |2] may be forbidden. The updating 
of the partial generating function depends most crucially on the state of vertex b and 
to a somewhat lesser extent on the states of the vertices a and d. The basic limitation 
on the allowed outputs are that conditions (2)-(4) must be enforced. In the following 
we shall briefly describe how the updating rules are derived. 

State of vertex b is 0. Since vertex b is empty all the outputs in figure El are possible. 
In the first output we insert no bonds. This is always allowed and no changes are 
made to the configuration. 

In the next three outputs we insert a single bond. This makes vertex b of degree 
one and is thus only allowed if there is at most one free end in the source. There are 
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further restrictions on the insertion of a bond to vertices a or d. Firstly if a vertex 
is touched (in state 3) we cannot insert a bond since this would result in a vertex 
of degree 3. Secondly if the vertex is a free end (in state 4) we join two free ends. 
This leads to the formation of a completed sub-graph and is only permitted if the 
resulting graph is a valid SAW. So the configuration cannot contain other pieces of 
a SAW and the only permissible states of other vertices in the intersection are 
and 3. If a valid SAW is created we multiply the source generating function by u 
(representing the new bond) before adding it to the total for the SAW generating 
function. 

In the last three outputs we insert a partial loop. Again there are restrictions on 
the insertion of bonds to vertices a and d. As before we cannot insert a bond to 
a vertex in state 3. Otherwise the first two outputs are always allowed. The last 
output is a little more complicated. If both vertices a and d are in state 4 we join 
two free ends and as before we check if the result is a valid SAW and if so add this 
partial generating function the SAW generating function (this time we multiply the 
source generating function by u 2 ). If vertex a is in state 1 and vertex d in state 2 
we cannot join the two vertices since this would result in a closed loop. 
After the insertion of new bonds we have to assign a state to vertex c and quite 
possibly change the states of vertices a and d (and perhaps the states of some other 
vertices in the target configuration). The state of vertex c will be (no bond), 
1 (lower loop-end), 2 (upper loop-end) or 4 (free end). Next we consider what 
happens to vertices a and d. When these vertices are empty in the source they can 
take the values just listed above in the target. If they are occupied in the source 
they either retain their state in the target (no bonds inserted) or change to state 3 
(a bond is inserted). In the latter case we may have to change the state of other 
vertices in the target. If we insert a free end and it joins a lower (upper) loop-end 
we must change the matching upper (lower) loop-end to a free end. Otherwise we 
may join two lower (upper) loop-ends and then we must change the matching upper 
(lower) loop-end of the inner most loop to the lower (upper) loop-end of the new 
joined loop. 

State of vertex b is 1. A lower end of a loop enters vertex b. If we insert no further 
bonds a new degree 1 vertex is created. As above this is only allowed provided 
the source has at most one free end. The matching upper loop-end becomes a free 
end. Otherwise the lower end has to be continued by inserting a single bond (partial 
loops cannot be inserted since this would make vertex b of degree 3) either to vertex 
c which becomes a state 1 vertex; to vertex a if not in state 3 or state 2 (closed loop 
would be formed); or to vertex d if not in state 3. Again we have to change the 
states of vertices a and d when a bond is inserted on these vertices. If the source 
state of the vertices was the target state becomes 1, otherwise the target state 
becomes 3 and as above we may need to change the state of other vertices as well. 

State of vertex b is 2. An upper end of a loop enters vertex b. If we terminate the 
loop-end a new degree 1 vertex is created. Again this is only allowed provided the 



Self-avoiding walks and polygons on the triangular lattice 



8 



source has at most one free end. The matching lower end of the loop becomes a 
free end. The upper end can always be continued to vertex c; to vertex d if it is 
not in state 3; and to vertex a provided it is not in state 3 or 1 (this would result 
in a closed loop). The state of the target vertices are changed as described above. 

State of vertex b is 3. This is the simplest situation. Vertex b is of degree 2 so no 
bonds can be inserted and only the output with all empty edges is allowed. The 
state of vertex c is and the states of all other vertices are unchanged. 

State of vertex b is 4. A free end is entering vertex b. If we insert no further bonds a 
partial walk is terminated at the vertex. This is only allowed if the resulting graph 
is a valid SAW and the source generating function is added to the SAW generating 
function. The free end can always be continued to vertex c and to vertices a and d 
if they are not in state 3. As before, if we join two free ends we check if it is a valid 
SAW and then add the partial generating function (multiplied by u) to the SAW 
generating function. Otherwise the target configuration is updated as described 
previously. 

SAP updating rules. The updating rules used when enumerating SAPs are essentially 
just a subset of the SAW rules. Obviously there are no degree 1 vertices in the SAP 
case so we can't insert a single bond if vertex b is empty. Likewise if vertex b is occupied 
we must continue the loop-end. Completed SAPs are formed by closing a loop (if there 
are no other loop-ends in the source). This happens when the local source configuration 
{abd} is {120} and we insert a bond from b to a, {102} and we insert a partial loop 
from a through b to d, or {012} and we insert a bond from b to d. 

2.1.2. Pruning The use of pruning to obtain more efficient TM algorithms was used for 
square lattice SAPs in [21] . Numerically it was found that the computational complexity 
was close to 2 n//4 , much better than the 3 n / 4 of the original approach. We have used 
similar techniques for the enumerations carried out for this paper. Pruning allows us 
to discard most of the possible configurations for large W because they only contribute 
at lengths greater than A" max , where iV max is the maximal length to which we choose to 
carry out our calculations. The value of iV max is limited by the available computational 
resources, be they CPU time or physical memory. Briefly pruning works as follows. 
Firstly, for each configuration we keep track of the current minimum number of steps 
N mT already inserted to the left of the boundary line in order to build up that particular 
configuration. Secondly, we calculate the minimum number of additional steps A" a dd 
required to produce a valid SAP or SAW. There are three contributions, namely the 
number of steps required to connect the loops and free ends, the number of steps needed 
(if any) to ensure that the SAW touches both the lower and upper border, and finally the 
number of steps needed (if any) to extend at least W edges in the length-wise direction 
(remember we only need rectangles with L > W). If the sum N cm + N a dd > N max 
we can discard the partial generating function for that configuration, and of course the 
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configuration itself, because it won't make a contribution to the count up to the lengths 
we are trying to obtain. 

There are no principal differences between pruning SAWs and SAPs though 
the detailed implementation is more complicated for the SAW case. We found it 
necessary to explicitly write subroutines to handle the three distinct cases where the 
TM configuration contains zero, one and two free ends, respectively. But in all cases 
we essentially have to go through all the possible ways of completing a SAW in order to 
find the minimum number of steps required. This is a fairly straightforward task though 
quite time consuming. 

2.1.3. Computational complexity The time T(n) required to obtain the number of 
polygons or walks of length n grows exponentially with n, T(n) oc A™. For the 
algorithm without pruning the complexity can be calculated exactly. Time (and 
memory) requirements are basically proportional to a polynomial (in n) times the 
maximal number of configurations, iVconf; generated during a calculation. When the 
boundary line is straight we can find the exact number of configurations. First look at 
the situation for SAPs when there are no free ends. The configurations correspond to 
2-coloured Motzkin paths ,9j], since we can map empty and touched vertices to the two 
colours of horizontal steps, vertices in the lower state to a north-east step, and vertices 
in the upper state to a south-east step. The number of such paths M n with n steps is 
easily derived from the generating function [0] 



which means that M n ~ 4" as n — > oo. Note that M n slightly over counts iVconf since 
configurations without a loop aren't permitted, but since there are only 2 W of these 
there is no change in the asymptotic growth. When the boundary line has a kink (such 
as in figure |TJ) iV Con f is no longer given exactly by M w — 2 W . However, it is obvious 
that N$(W + 1) < iVconf < Ns(W) so we see that asymptotically iVconf grows like 4 W . 
Since a calculation using rectangles of widths < W yields the number of SAPs up to 
n = 2W + 1 it follows that the complexity of the algorithm is T(n) oc A n with A = 2. 

The number of transfer matrix configurations in the unpruned SAW algorithm is 
simply obtained by inserting 0, 1 or 2 free ends into a 2-coloured Motzkin path and 
eliminating the paths corresponding to a configurations with only empty or touched 
vertices. In this case a calculation using rectangles of widths < W yields the number of 
SAWs up to n = W it follows that the complexity of the algorithm is T{n) oc A n with 



The pruned algorithm is much too difficult to analyse exactly. So all we can 
give is a numerical estimate of the growth in the number of configurations with n. 
That is obtained by just running the algorithm and measuring the maximal number of 
configurations generated for different values of n. From the actual numbers it appears 
that for the SAP case increasing n by 2 increases the number of configurations by slightly 
less than a factor of 2. This would mean that for the pruned SAP algorithm A p ~ a/2. 




(2.2) 



■n 



A = 4. 
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In the SAW case it appears that increasing n by 4 increases the number of configuration 
by a factor close to 5. So for the pruned SAW algorithm X p w \fh = 1.495 .... So once 
again pruning results in much more efficient algorithms. 

2.1.4- Parallelisation The transfer-matrix algorithms used in the calculations of the 
finite lattice contributions are eminently suited for parallel computation. The bulk of the 
calculations for this paper were performed on the facilities of the Australian Partnership 
for Advanced Computing (APAC). The APAC facility is an HP Alpha server cluster 
with 125 ES45's each with four 1 Ghz chips for a total of 500 processors in the compute 
partition. Each server node has at least 4 Gb of memory. Nodes are interconnected by 
a low latency high bandwidth Quadrics network. 

The most basic concern in any efficient parallel algorithm is to minimise the 
communication between processors and ensure that each processor does the same 
amount of work and uses the same amount of memory. In practice one naturally has to 
strike some compromise and accept a certain degree of variation across the processors. 

One of the main ways of achieving a good parallel algorithm using data 
decomposition is to try to find an invariant under the operation of the updating rules. 
That is we seek to find some property of the configurations along the boundary line which 
does not alter in a single iteration. The algorithm for the enumeration of SAPs and 
SAWs are quite complicated since not all possible configurations occur due to pruning 
and an update at a given set of vertices might change the state of a vertex far removed, 
e.g., when two lower loop-ends are joined we have to relabel one of the associated 
upper loop-ends as a lower loop-end in the new configuration. However, there is still 
an invariant since any vertex not directly involved in the update cannot change from 
being empty to being occupied and vice versa, likewise a touched vertex will remain 
unchanged. That is only the kink vertices can change their occupation or touched status. 
This invariant allows us to parallelise the algorithm in such a way that we can do the 
calculation completely independently on each processor with just two redistributions of 
the data set each time an extra column is added to the lattice. We have already used 
this scheme for SAPs jTHj and lattice animals [TH] and refer the interested reader to 
these publications for further detail. Our parallelisation scheme is also very similar to 
that used by Conway and Guttmann |H1 IT3]. 

2.2. Metric properties and area-weighted moments 

The generalisation of the algorithm required to calculate metric properties and area- 
weighted moments has been described in detail in [171120] in the square lattice case. Only 
some minor adjustments are needed in order to apply these ideas to metric properties 
on the triangular lattice (no changes are needed for the area- weighted moments). We 
have represented the triangular lattice as a square lattice with extra edges along one of 
the main diagonals in a unit cell. A point (s,t) on the square lattice is the point (x,y) 
on the triangular lattice where x = s + \t and y = ^t. As shown in [T7J |20] calculation 
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of metric properties involves summation over products of the x and y coordinates of 
the distance vectors. To be explicit we define the radius of gyration according to the 
vertices of the SAW. Note that the number of vertices is one more than the number of 
steps. The radius of gyration of n + 1 points at positions is 

(n + l) 2 (R 2 g ) n = - r,) 2 = n J^i^ + Vi) ~ 2 J^&Xj + y iVj ). (2.3) 

i>j i i>j 

From the triangular lattice coordinates we see that both XiXj and y^yj carry a factor 
| so in order to ensure that we get integer coefficients we multiply by 4 and the algorithm 
will thus calculate the coefficients 4(n+ l) 2 c n (R 2 g ) n . In order to do this we maintain five 
partial generating functions for each possible boundary configuration, namely 

• C(u), the number of (partially completed) SAWs. 

• X 2 (u), the sum over SAWs of the squared components of the distance vectors. 

• X g (u), the sum of the x-component of the distance vectors. 

• Y g (u), the sum of the y-component of the distance vectors. 

• XY g (u), the sum of the 'cross' product of the components of the distance vectors, 
that is, Y^i>j( x i x j +ViVj)- 

As the boundary line is moved to a new position each target configuration S might 
be generated from several sources S' in the previous boundary position. The partial 
generating functions are updated as follows, with (s, t) being the coordinates of the 
newly added vertex on the square lattice: 



C(u,S) = ^u n 'C{u,S'), 

S' 

x >> s ) = J2 un ^ x >^ s '^ +6 ^+tr+3t 2 )c(u,s')], 

S' 

X g (u,S) =J2 un '{ X 9 ( u ^) + S g (2s + t)C(u,S% (2.4) 

S' 

Y g (u,S) = Y,u n '[Y g (u,S) + 5gtC(u,S% 

S' 

XY g {u, S) = ^/[iy 9 (n, S') + 5 g (2s + t)X g (u, S') + 5g3tY g (u, S')}, 

S' 

where n' is the number of steps added to the partial SAW. 8 g — if the new vertex is 
empty (has degree 0) and 8 g — 1 if the new vertex is occupied (has degree > 0). 
The updating rules for the other metric properties are generalised similarly. 



2.3. Enumeration results 

We calculated the number of polygons up to perimeter 60, while the radius of gyration 
and first 10 area-weighted moments were obtained up to perimeter 58. We calculated 
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the number of SAWs, their mean-square radius of gyration, mean-square end-to-end 
distance, and the mean-square distance of monomers from the end points. These 
quantities were obtained for walks up to length 40. The calculations required up to 
35Gb of memory using up to 40 processors at a time and in total we used about 15000 
CPU hours. 

In tabled we list the number of SAPs and their radius of gyration while in table [21 
we list the series for the SAW problem. These series and those for the area-weighted 
moments are available at www.ms.unimelb.edu.au/~iwan. 

3. Analysis of the series 

The series studied in this paper have coefficients which grow exponentially, with sub- 
dominant term given by a critical exponent. The generic behaviour is G(u) = Yl n 9nU n ~ 
(1 — u/u c )~^, and hence the coefficients of the generating function g n ~ //Vi^ -1 , where 
H — l/u c . To obtain the singularity structure of the generating functions we used the 
numerical method of differential approximants [12J. Our main objective is to obtain 
accurate estimates for the connective constant \x and confirm numerically the exact 
values for the critical exponents a, 7 and v. Estimates of the critical point and critical 
exponents were obtained by averaging values obtained from second and third order 
inhomogeneous differential approximants. The error quoted for these estimates reflects 
the spread (basically one standard deviation) among the approximants. Note that these 
error bounds should not be viewed as a measure of the true error as they cannot include 
possible systematic sources of error. 

Once the exact values of the exponents have been confirmed we turn our attention 
to the "fine structure" of the asymptotic form of the coefficients. In particular we are 
interested in obtaining accurate estimates for the amplitudes. We do this by fitting 
the coefficients to the form given by (jl.lajl - JT.lgj ). In this case our main aim is to 
test the validity of the predictions for the amplitude combinations mentioned in the 
Introduction. 

3.1. Self-avoiding polygons 

The expected behaviour of the mean-square radius of gyration Jl.l/D and moments of 
area fll.lg| ) of SAPs results in the following predictions for the generating functions 
which we study: 

00 00 

n]{u) = J2^Pn(R 2 )nU n = Ys T nU n OC (1 - U^ M \ (3.1) 

n=3 n=3 

oo oo 

V {k \u) = Y,Pn(a k ) n u n = J]a<V oc (1 - Uf i) 2 ^ a+2kv \ (3.2) 

n=3 n=3 

where we have taken into account that the smallest polygon has perimeter 3. Thus we 
expect these series to have a critical point at u = u c = 1/fi, and as stated previously 





Table 1. 


The number, p n , of embeddings of n-step polygons on the triangular lattice and their radius of gyration. 


n 


Pn 


p n n z {R z ) n 


n 


Pn 


p n n z {R z ) n 


3 


2 


6 


32 


2692047018699717 


25886228326621869696 


4 


3 


24 


33 


10352576717684506 


110846359749047031012 


5 


6 


102 


34 


39902392511347329 


474213717578995665624 


6 


15 


468 


35 


154126451419554156 


2026979522666735966994 


7 


42 


2172 


36 


596528356905096920 


8657009828812246231296 


8 


123 


9978 


37 


2313198287784319026 


36944420238568755696168 


9 


380 


45816 


38 


8986249863419780682 


157546885404468362432148 


10 


1212 


208686 


39 


34969337454759091232 


671378005865890422968520 


11 


3966 


944766 


40 


136301962040079085257 


2859142640844460643187642 


12 


13265 


4253484 


41 


532093404471021533628 


12168301979788445465498400 


13 


45144 


19046580 


42 


2080235431107538787148 


51756227545091330753357904 


14 


155955 


84891654 


43 


8144154378525048003270 


220011744770726296282498056 


15 


545690 


376756392 


44 


31927176350778729318192 


934740492588407244896782986 


16 


1930635 


1665684774 


45 


125322778845662829008494 


3969252848247139670605665948 


17 


6897210 


7338822888 


46 


492527188641409773340797 


16846468953704095289170900908 


18 


24852576 


32233105398 


47 


1937931188484341585677962 


71466199766730550647612342396 


19 


90237582 


141171369444 


48 


7633665703654150673637363 


303035054640652779166447899354 


20 


329896569 


616694403366 


49 


30101946001283232799847562 


1284380183482800747257353493532 


21 


1213528736 


2687630355198 


50 


118823919397444557546535851 


5441398704214816650431847144246 


22 


4489041219 


11687756315940 


51 


469508402822449711313115200 


23043633507948438933442640818176 


23 


16690581534 


50726031551790 


52 


1856933773092076293566747007 


97548735673726189271333029096494 


24 


62346895571 


219753786787212 


53 


7351015093472721439659392448 


412789876403022674873495520537906 


25 


233893503330 


950403133411176 


54 


29126027071450640626653986531 


1746140617537848477455116275581178 


26 


880918093866 


4103923685277414 
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27 
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56 
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28 


12630175810968 


76195720234557276 


57 
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29 


48056019569718 
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58 


7237327695683743010999188700157 


557209110506518251250962658184410206 


30 


183383553173255 


1407546930663067986 


59 


28789332223533619621001538109842 




31 


701719913717994 


6039368800117995984 


60 


114602547490254934327469368968190 





Table 2. The number, c n , of embeddings of n-step self-avoiding walks on the triangular lattice and their radius of gyration, end-to-end 
distance and distance of monomers from the end-points. 



Tl 






l( n + l) 2 Cn(R 2 ) n 


l( n + l)c n (R 2 )„ 


1 


6 


l 


1 


1 


2 


30 


12 


22 


17 


3 


138 


97 


282 


178 


4 


618 


654 


2778 


1476 


5 


2730 


3977 


23305 


10667 


6 


11946 


22624 


175194 


70359 


7 


51882 


122821 


1215740 


434708 


8 


224130 


644082 


7939156 


2557166 


9 


964134 


3288739 


49422491 


14477823 


10 


4133166 


16440648 


295993366 


79492861 


11 


17668938 


80783857 


1717056604 


425633898 


12 


75355206 


391310240 


9697408184 


2231674940 


13 


320734686 


1872763387 


53533130211 


11494836257 


14 


1362791250 


8870963422 


289769871988 


58310378811 


15 


5781765582 


41647686501 


1541876281342 


291901836462 


16 


24497330322 


194014270964 


8081886977224 


1444405248178 


17 


103673967882 


897639074623 


41801262603145 


7074419785415 


18 


438296739594 


4127904278590 


213650877117460 


34334678700977 


19 


1851231376374 


18879838654237 


1080407596025856 


165283451747722 


20 


7812439620678 


85930246593928 


5411153165106856 


789827267540498 


21 


32944292555934 


389382874004291 


26865804448156781 


3749241090582031 


22 


138825972053046 


1757383045067340 


132328831054383256 


17689855417349797 


23 


584633909268402 


7902553525660965 


647064413113509344 


83004601828121876 


24 


2460608873366142 


35417121500633314 


3142945284616515512 


387503899136724032 


25 


10350620543447034 


158241760294727837 


15172247917136636793 


1800616777561080887 


26 


43518414461742966 


705008848574456242 


72826367061554681960 


8330920471773661365 


27 


182885110185537558 


3132749279518281223 


347722481262776946768 


38390978707292879316 


28 


768238944740191374 


13886614514918779812 


1652126117509776447678 


176259763248055992656 


29 


3225816257263972170 


61415827107198652263 


7813839241496101017943 


806446563482615080995 


30 


13540031558144097474 


271046328280157919578 


36798230598686798952874 


3677867046530479086571 


31 


56812878384768195282 


1 193838903060544883615 


172603075240086498030932 


16722626138383080469074 


32 


238303459915216614558 


5248569464050058190772 


806559315077883801952302 


75819788411079420147060 


33 


999260857527692075370 


23034474248167644819305 


3755672941408238341746325 


342850281196290726391195 


34 


4188901721505679738374 


100925879660029490332616 


17429779928912903943728776 


1546457563237807336247617 


35 


17555021735786491637790 


441524252843364233569911 


80636231608943399450377104 


6958970268567678359172166 


36 


73551075748132902085986 


1928731794198995523104424 


371943975622752362856339418 


31245121332848941331142166 


37 


308084020607224317094182 


8413734243045682304542891 


1710813401690158618688146075 


139991577634597301110308061 


38 


1290171266649477440877690 


36655327788277288494374240 


7848181414990001769700643892 


625968026891459936611240307 


39 


5401678666643658402327390 


159494618902280757690831541 


35911648943670829119431170002 


2793684462154188994667777314 


40 


22610911672575426510653226 


693174559672551318610401776 


163929038497681452701025717812 


12445679176337664122926617782 
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Table 3. Estimates for the critical point u c and exponent 2 — a obtained from second 
and third order differential approximants to the triangular lattice SAP generating 
function. L is the order of the inhomogeneous polynomial. 



L 


Second order DA 


Third order DA 




u c 2 — a 


u c 2 — a 



0.2409175671(28) 

2 0.2409175709(14) 

4 0.2409175714(27) 

6 0.2409175707(29) 

8 0.2409175724(44) 

10 0.2409175717(39) 



1.5000142(45) 

1.5000076(29) 

1.5000061(56) 

1.5000075(58) 

1.500003(10) 

1.5000051(83) 



0.2409175706(62) 
0.2409175716(30) 
0.2409175699(40) 
0.2409175712(29) 
0.2409175662(80) 
0.2409175704(22) 



1.500006(12) 

1.5000071(58) 

1.5000078(63) 

1.5000065(57) 

1.500012(14) 

1.5000083(41) 



the exponents a = 1/2 and v = 3/4. 

3.1.1. The SAP generating function In Table 0] we have listed the results from our 
series analysis of the SAP generating function. It is evident that the estimates for the 
critical exponent is in complete agreement with the expected value 2 — a = 3/2. Based 
on the estimates we find that u c = 0.24091757(1). We found no evidence that the SAP 
generating function had any other singularities. 

If we take the conjecture a = 1/2 to be true we can obtain a refined estimate for 
the critical point u c . In figure |H] we have plotted estimates for the critical exponent 
2 — a against the number of terms used by the approximant and against estimates 
for the critical point u c , respectively. Each dot represents estimates obtained from a 
third order inhomogeneous differential approximant. The order of the inhomogeneous 
polynomial was varied from to 10. As can be seen from the left panel the estimates 
for the critical exponent clearly include the exact value and appear to settle down as 
the number of terms increases (though there is a hint of a downwards trend). From 
the right panel we observe that the estimates cross the line 2 — a = 3/2 at a value 
u c ~ 0.2409175745. Based on this analysis we adopt the value u c = 0.2409175745(15) 
and thus // = 4.150797226(26) as our final estimates. 

3.1.2. The radius of gyration Table 0] lists the results of our series analysis for 
the SAP radius of gyration generating function. It is evident that the estimates 
for the critical exponent as obtained from third order differential approximants is in 
complete agreement with the expected behaviour. The estimates from the second 
order approximants are generally slightly lower than the expected value. One would 
generally expect third order differential approximants to be superior since they are better 
suited to represent complicated functional behaviour. We take this as clear numerical 
confirmation that v = 3/4. 

3.1.3. Area-weighted moments In passing we only briefly mention that our analysis of 
the generating functions V^ k \u) for area- weighted SAPs with k < 10 clearly confirmed 
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Figure 3. Plots of estimates from third order differential approximants for 2 — a vs., 
respectively, the number of terms used by the differential approximants (left panel) 
and u c (right panel). 



Table 4. Estimates for the critical point u c and exponent a + 2v of the SAP radius 
of gyration generating function. 



L Second order DA Third order DA 





u c 


a + 2v 


u c 


a + 2v 





0.24091726(11) 


1.99885(27) 


0.24091761(30) 


2.0001(14) 


2 


0.24091727(14) 


1.99891(26) 


0.24091728(33) 


1.99925(64) 


4 


0.240917246(95) 


1.99881(14) 


0.24091713(30) 


1.99889(36) 


G 


0.240917269(87) 


1.99884(14) 


0.24091741(19) 


1.99935(65) 


8 


0.240917239(73) 


1.99879(11) 


0.24091743(24) 


1.99947(78) 


10 


0.240917281(96) 


1.99888(16) 


0.24091737(25) 


1.99932(70) 



the expected values, 2 — (a + 1kv\ for the critical exponents. Suffice to say that the 
estimates range from 0.0005(8) for k — 1 to —13.503(2) for for k = 10. 

3.2. Self-avoiding walks 

From the expected behaviour (jl.lajl of c n and the metric properties of SAWs (jl.lcjl - 
(ll.le|) we get that the generating functions: 

oo 

C(u) = ^c n u n oc (1 -M/i)~ 7 , (3.3) 

n=l 

oo 

TII{U) = Y, C n(Rl)nU n (X (1 -U/l)-^\ (3.4) 
n=l 
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Table 5. Estimates for the critical point u c and exponent 7 obtained from second and 
third order differential approximants to the square lattice SAW generating function. 



L 


Second order DA 


Third order DA 




u c 7 


u c 7 



0.240917491(34) 

2 0.240917529(37) 

4 0.240917529(42) 

6 0.240917524(27) 

8 0.240917523(28) 

10 0.240917513(31) 



1.343637(42) 
1.343677(36) 
1.343682(47) 
1.343673(27) 
1.343668(35) 
1.343662(29) 



0.240917538(21) 
0.240917537(13) 
0.240917534(30) 
0.240917545(24) 
0.240917543(23) 
0.240917530(22) 



1.343687(23) 
1.343686(22) 
1.343682(32) 
1.343693(25) 
1.343692(27) 
1.343679(25) 



K( U ) = + l) 2 Cn(^)n«" OC (1 - W )-^ +2 \, 

n=l 

00 

KI(U) = J> + VCniRDnU" (X (1 - U^ + ^ +1 \ 
n=l 

where the exponents 7 = 43/32 and v — 3/4. 



(3.5) 
(3.6) 



3.2.1. The SAW generating function In Table El we list estimates of the critical point 
u c and exponent 7 from the series for the triangular lattice SAW generating function. 
The estimates were obtained by averaging over those approximants which used at 
least the first 32 terms of the series. Based on these estimates we conclude that 
u c = 0.24091753(8) and 7 = 1.34368(6). The estimate for u c is compatible with the 
more accurate estimate u c = 0.2409175745(15) obtained above from the analysis of the 
SAP generating function. The analysis adds further support to the already convincing 
evidence that the critical exponent 7 = 43/32 = 1.34375 exactly. However, we do 
observe that both the central estimates for both u c and 7 are systematically slightly 
lower than the expected values. 

As for the SAP case we find it useful to plot the behaviour of the estimates for the 
critical exponent 7 against both u c and the number of terms used by the differential 
approximants. This is done in figure El Each dot represents estimates obtained from a 
second order inhomogeneous differential approximant. From the left panel we observe 
that the estimates of 7 exhibits a certain upwards drift as the number of terms increases. 
So the estimates have not yet settled at their limiting value, but there can be no doubt 
that the predicted exact value of 7 is fully consistent with the estimates. In the left 
panel we observe that the (u c , 7)-estimates fall in a narrow range. Note that this range 
does not include the intersection point between the exact 7 and the precise u c estimate. 
This is probably a further reflection of the lack of 'convergence' to the true limiting 
values. 
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Figure 4. Plots of estimates from second order differential approximants for 7 vs., 
respectively, the number of terms used by the differential approximants (left panel) 
and u c (right panel). 



Table 6. Estimates for the critical point u c and critical exponents as obtained 
from third order differential approximants to the generating functions for the metric 
properties of SAWs. 



L 


Kc{ 


A 


K g (u) 






U c 


■y + 2u 


U c 


7 + 2^ + 2 


U c 


7 + 2J/ + 1 





0.240917330(86) 


2.84307(36) 


0.240917594(53) 


4.843619(70) 


0.240917324(92) 


3.84296(17) 


2 


0.240917298(62) 


2.84295(12) 


0.240917600(53) 


4.843626(72) 


0.24091715(22) 


3.84270(35) 


4 


0.240917249(39) 


2.84272(35) 


0.240917605(62) 


4.843631(81) 


0.24091722(23) 


3.84281(41) 


6 


0.240917311(71) 


2.84295(16) 


0.240917578(71) 


4.843590(99) 


0.24091732(17) 


3.84299(34) 


8 


0.240917328(52) 


2.842938(73) 


0.240917616(67) 


4.843646(89) 


0.240917304(43) 


3.842922(80) 


10 


0.240917373(99) 


2.84303(19) 


0.240917612(57) 


4.84365(10) 


0.240917276(98) 


3.84285(18) 



3.2.2. The metric properties Finally, we turn our attention to the metric properties of 
SAWs. In Table El we have listed the estimates for u c and critical exponents obtained by 
an analysis of the associated generating functions (|H.4j) - (|H.fi)L The estimates from the 
radius of gyration u c = 0.2409176(1) and 7 + 2z/ + 2 = 4.84365(20) are in full agreement 
with the more accurate SAP estimate for u c and the predicted exact exponent value 
7 + 2u + 2 = 155/32 = 4.84375. The analysis of the generating functions for the end-to- 
end distance and monomer distance yield estimates of u c a little below the expected value 
and likewise the exponent estimates 2.8430(5) and 3.8429(5) are a somewhat smaller that 
the exact values 7 + 2u + 2 = 91/32 = 2.84375 and 7 + 2u + 2 = 123/32 = 3.84375, 
respectively. We are fully convinced that this is because the series are not long enough 
to allow the exponent estimates to settle at the true limiting values, as was also the case 
for the SAW generating function as shown in figure 0] 
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3.3. Amplitudes 

The asymptotic form of the coefficients p n of the generating function of square lattice 
SAPs has been studied in detail previously [HJ El ESI- There is now clear numerical 
evidence that the leading correction-to-scaling exponent for SAPs and SAWs is Ai = 
3/2, as predicted by Nienhuis (2E1I2Z1- As argued in [Hj this leading correction term 
combined with the 2 — a = 3/2 term of the SAP generating function produces an 
analytic background term. Indeed in the previous analysis of SAPs there was no 
sign of non-analytic corrections-to-scaling to the generating function (a strong indirect 
argument that the leading correction-to-scaling exponent must be half-integer valued). 
One therefore finds that asymptotically p n behaves as 



This form was confirmed with great accuracy in j2U CHj- Estimates for the leading 
amplitude B can thus be obtained by fitting p n to the form (|3.7[) using an increasing of 
number of correction terms. As in we find it useful to check the behaviour of the 
estimates by plotting the results for the leading amplitude vs. 1/n (see figure E}, where 
p n is the last term used in the fitting. In addition we also wish to check the sensitivity 
of the procedure to small changes in the value of \x. Clearly the amplitude estimates 
in top panels are quite well converged. Notice that as more correction terms are added 
the estimates exhibit less curvature and the slope becomes less steep. This is very 
strong evidence that ()3.7j) indeed is the correct asymptotic form of p n . The estimates 
shown in the bottom panels are not so well behaved. Those in the left panel are not 
monotonic and after initially decreasing they start to increase with n. The estimates in 
the right panel while monotonic have much steeper slopes and the slopes do not appear 
to change much as more correction terms are used. We think this is strong evidence 
that \i = 4.150797226 is very close to the true value. Based on the plots in the top right 
panel we estimate that B = 0.2639393(1). 

The asymptotic form of the coefficients r n in the generating function for the radius 
of gyration was studied in J7|. The generating function (|3.1|) has critical exponent 
— (a + 2u) = —2, so the leading correction-to-scaling term no longer becomes part of 
the analytic background term. We thus use the following asymptotic form: 



We find BF = 0.013710(1). This is in full agreement with the predicted exact value |I] 
= 3§^ aa o = = 0.013710424. . ., where, for the triangular lattice, a = 1 and 
a = V3/2. Combining the exact expression for BF with the more accurate estimate 
for B given above we estimate that F = 0.05194537(2). 




(3.7) 




(3.8) 
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Figure 5. Estimates of the leading amplitude B plotted against 1/n using different 
number of terms in the asymptotic expansion (|3.7I) . In the top panels we use the value 
H — 4.150797226. The right panel just gives a more detailed look at the data shown 
in the left panel. In the bottom panels we use two different values [i = 4.15079720 
(left panel) and /i = 4.15079725 (right panel) chosen to be at the extreme ends of our 
error-bounds for /x. 



The amplitudes of the area-weighted moments were studied in 
coefficients to the assumed form 



np n (a ) n ~ n n n 



n2ku+a.—2 



k\ 



ft 



i>0 



We fitted the 



(3-9) 



where the amplitude G k = G^ k 'B/k\ is related to the amplitude defined in equation 
( l.lg ). The scaling function prediction for the amplitudes G k is [30J 

b 2k (3k-2)\ ^ n2k b 2k +i 



Go k B 2k ~ 1 



4vr 3fc (6fc-3)!' + (3k)h 3k + 1 2 6k + 2 ' 

where the numbers b k are given by the quadratic recursion 

j fe-i 

b k + (3k - 4)6 fc _i + - ^6 fc _ r 6 r = 0, b = 1. 



(3.10) 



(3-11) 



r=l 
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Table 7. Predicted exact values for universal amplitude combinations and estimates 
from enumeration data for square, hexagonal and triangular lattice polygons. 



Amplitude 


Exact value 




Square 




Hexagonal 


Triangular 


B 


unknown 




0.56230130(2) 




1.27192995(10) 


0.2639393(2) 


Gi 


0.7957747 x 1CT 


1 


0.795773(2) x 10~ 


1 


0.795779(5) x 10" 1 


0.795765(10) x 10' 1 


G 2 B 


0.3359535 x 10' 


2 


0.335952(2) x 10~ 


2 


0.335957(6) x 10~ 2 


0.335947(5) x 10~ 2 


G3B 2 


0.1002537 x 10~ 


3 


0.100253(1) x 10~ 


3 


0.100255(3) x 10~ 3 


0.100251(4) x 10~ 3 


G4.B 3 


0.2375534 x 10" 


5 


0.237552(2) x 10" 


5 


0.237557(7) x 10~ 5 


0.237547(6) x 10~ 5 


G5B 4 


0.4757383 x 10" 


7 


0.475736(3) x 10~ 


7 


0.475749(10) x 10" 7 


0.475724(15) x 10" 7 


G 6 B 5 


0.8366302 x 10" 


9 


0.836624(5) x 10" 


9 


0.836652(10) x 10" 9 


0.83660(2) x 10~ 9 


G7B 6 


0.1325148 x 10" 


10 


0.132514(2) x 10" 


10 


0.132519(5) x 10" 10 


0.132511(5) x 10" 10 


GgB 7 


0.1924196 x 10" 


12 


0.192418(2) x 10- 


12 


0.192426(8) x 10~ 12 


0.192419(8) x 10" 12 


G9B 8 


0.2594656 x 10" 


14 


0.259464(2) x 10" 


14 


0.259472(12) x 10" 14 


0.25948(4) x 10~ 14 


G W B 9 


0.3280633 x 10" 


16 


0.328062(4) x 10~ 


l(i 


0.328051(15) x 10" 16 


0.32812(5) x lO" 16 



We obtained [3T] the results for the amplitude combinations listed in Tabled It is clear 
that the estimates for the first 10 area weighted moments are in excellent agreement 
with the predicted exact values. 

The amplitude ratios D/C and E/C were estimated by direct extrapolation of the 
relevant quotient sequence, using the following method |2H1 : Given a sequence {a n } 
defined for n > 1, assumed to converge to a limit a m with corrections of the form 



[l+b/n+. . .), we first construct a new sequence {h n } defined by h n = Yi 



m=l ^m- 



The generating function h n x n ~ (1 — a^x) Estimates for and the parameter 

b can then be obtained from differential approximants. In this way, we obtained the 
estimates [I] D/C = 0.140296(6) and E/C = 0.439649(9). These amplitude estimates 
leads to a high precision confirmation of the CSCPS relation (jl.2J) H = 0.000036(34). 

The amplitudes of the SAW generating function and the metric properties were 
also studied in pQ by fitting of the coefficients to the assumed form 



A + J2^/n 1+l/2 
c n (R 2 e ) n /6~n n n 



i>0 



(n + l) 2 c n (^) n /6~/i n n 



AC/6 + J2 a i/ nl+l/2 

i>0 

AD/Q + J2 a i/ nl+l/2 



n 7+2Z/+1 



i>0 



(n + l)c n «) n /6~/i n n 



,"^7+2^ 



AE/6 + 



/2 



i>0 



(3.12) 
(3.13) 
(3.14) 
(3.15) 



In figure |H1 we have plotted the estimates for the leading amplitudes against 1/n 
while varying the number of terms used in fitting to the expressions given above. From 
these we estimate that A = 1.183966(1), AC/6 = 0.140380(5), AD/6 = 0.019696(3), 
and AE/6 = 0.061715(5). The estimate for A is the same as that obtained in £Q 
while the remaining amplitude estimates are a little lower and have smaller error-bars 
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Figure 6. Estimates of the leading amplitude A for SAWs, AC/6 for the end-to-end 
distance, AD/ 6 for the radius of gyration, and AE/6 for the monomer distance from 
the end-points, plotted against 1/n using varying number of terms in the asymptotic 
expansion. 

that those quoted in pQ. The main reason is that here we are only interested in the 
leading amplitudes and base our estimates on the fits using 6-8 terms, while in [T] 
estimates for sub-leading amplitudes were also considered and required to be stable and 
consequently only fits with up to 4 terms were considered. For the metric amplitudes 
we thus obtain the estimates C = 0.71140(3), D = 0.099814(15), and E = 0.31275(3). 
For the universal amplitude ratios we get D/C = 0.14030(2) and E/C = 0.43963(5). 
We note that these estimates of the amplitude ratios are fully consistent with the more 
accurate estimates given above. This gives us further confidence that this method 
for obtaining amplitude estimates is valid. In particular, it appears, that in order to 
estimate the leading amplitude, we do not have to insist that estimates for sub-leading 
amplitudes be well converged. The smaller error-bars obtained from the fits using 6-8 
terms thus appear soundly based. Naturally, some readers might wish to take a more 
cautious approach. 
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Table 8. Estimates of universal amplitude combinations on some two-dimensional 
lattices. 



Lattice 


D/C 


E/C 


BC/aao 


H 


Square PEH] 


0.140299(6) 


0.439647(6) 


0.216835(15) 


-0.000024(28) 


Triangular pQ 


0.140296(6) 


0.439649(9) 


0.216823(10) 


-0.000036(34) 


Honeycomb |22I 


0.1403(1) 


0.4397(2) 


0.2170(3) 


-0.00013(67) 


Kagome [53 


0.140(1) 


0.440(1) 


0.2144(25) 


-0.0015(47) 



In Table [H] we have summarised estimates of various universal amplitude 
combinations, obtained from the work reported in this paper and elsewhere. As can 
be seen the estimates for all lattices are in perfect agreement strongly confirming the 
universality of the various combinations. 

4. Summary and conclusion 

We have presented both improved and parallel algorithms for the enumeration of self- 
avoiding polygons and walks on the triangular lattice. These algorithms have enabled us 
to obtain polygons up to perimeter length 60, their radius of gyration and area- weighted 
moments up to perimeter 58, while for self-avoiding walks to length 40 we calculated the 
number of walks as well as the metric properties of mean-square end-to-end distance, 
mean-square radius of gyration and the mean-square distance of a monomer from the 
end points. 

The analysis of the polygon series enabled us to obtain a very precise estimate 
for the connective constant fi = 4.150797226(26). We confirmed to a very high degree 
of accuracy the predicted exponent values a = 1/2, 7 = 43/32 and v = 3/4. We 
noticed that, as is the case for the square lattice problem, the SAW asymptotics is worse 
behaved than the SAP asymptotics, i.e., estimates for /1 and the critical exponents are 
at least an order of magnitude more accurate in the SAP case. It quite is possible 
that this behaviour is due to the leading correction-to-scaling exponent A = 3/2. In 
the SAP case this correction simply becomes part of the analytic background term and 
the SAP generating function is therefore simpler since it only has analytic corrections 
to scaling. We also obtained very accurate estimates for the leading amplitude of the 
sequence p n of SAP coefficients B = 0.2639393(1) and using the exact expression for the 
amplitude combination BF we find F = 0.05194537(2). Our data for the area-weighted 
moments was used [SI] to confirm the correctness of theoretical predictions for the 
values of the amplitude combinations G^B k ~ 1 . Finally we obtained accurate estimates 
for the critical amplitudes A = 1.183966(1), C = 0.71140(3), D = 0.099814(15), and 
E = 0.31275(3). The estimate for the ratio C/F = 13.6952(5) is in very good agreement 
with the theoretical estimate C/F rs 13.70 [Hj. The amplitude estimates led to a high 
precision confirmation of the CSCPS relation ()1.2|) H = 0. 
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E-mail or WWW retrieval of series 

The series for the problems studied in this paper can be obtained via e-mail by sending 
a request to I.Jensen@ms.unimelb.edu.au or via the world wide web on the URL 
http://www.ms. unimelb.edu.au/~iwan7] by following the relevant links. 
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